Superconductivity in CeCa explained. 



o 
o 

(N 



o 
o 

I 

Oh 



I 

O 

o 



(N 
> 

00 

o 

O 
in 
o 



I 

O 

o 



Matteo Calandra and Francesco Mauri 
Institut de Mineralogie et de Physique des Milieux condenses, 
case 115, 4 place Jussieu, 75252, Paris cedex 05, France 
(Dated: February 2, 2008) 

Using density functional theory we demonstrate that superconductivity in CeCa is due to 
a phonon-mediated mechanism with electron-phonon coupling A = 0.83 and phonon-frequency 
logarithmic-average (cj) = 24.7 meV. The calculated isotope exponents are Q(Ca) — 0.24 and 
a(C) = 0.26. Superconductivity is mostly due C vibrations perpendicular and Ca vibrations paral- 
lel to the graphite layers. Since the electron-phonon couplings of these modes are activated by the 
presence of an intercalant Fermi surface, the occurrence of superconductivity in graphite intercalated 
compounds requires a non complete ionization of the intercalant. 

PACS numbers: 74.70.Ad, 74.25. Kc, 74.25. Jb, 71. 15. Mb 



Graphite intercalated compounds (GICs) were first 
synthesized in 1861 Q but only from the 30s a systematic 
study of these systems began. Nowadays a large number 
of reagents can be intercalated in graphite 100)0- 
Intercalation allows to change continuously the proper- 
ties of the pristine graphite system, as it is the case for the 
electrical conductivity. The low conductivity of graphite 
can be enhanced to obtain even larger conductivities than 
Copper [3[. Moreover at low temperatures, intercalation 
can stabilize a superconducting state The discovery 
of superconductivity in other intercalated structures like 
MgB2 y and in other forms of doped Carbon (diamond) 
has renewed interest in the field. 

The first discovered QIC superconductors were alkali- 
intercalated compounds 6] (CgA with A= K, Rb, Cs with 
Tc < 1 K). Synthesis under pressure has been used to 
obtain metastable GICs with larger concentration of al- 
kali metals (CgK, C3K, C4Na, C2Na) where the high- 
est Tc corresponds to the largest metal concentration, 
Tc(C2Na)=5 K 7]. Intercalation involving several stages 
have also been shown to be superconducting 

[113 (the 

highest Tc — 2.7 K in this class belongs to KTI1.5C4). In- 
tercalation with rare-earths has been tried, CgEu, CgCm 
and CgTm are not superconductors, while recently it has 
been shown that CeYb has a Tc = 6.5 K jlOlJ. Most 
surprising superconductivity on a non-bulk sample of 
CgCa was also discovered 10]. The report was confirmed 
by measurements on bulk CgCa poly-crystals 11] and a 
Tc ~ 11.5 K was clearly identified. At the moment CgYb 
and CgCa are the GICs with the highest Tc. It is worth- 
while to remember that elemental Yb and Ca are not 
superconductors. 

Many open questions remain concerning the origin of 
superconductivity in GICs. (i) All the aforementioned in- 
tercalants act as donors respect to graphite but there is 
no clear trend between the number of carriers transferred 
to the Graphene layers and TcQ- What determines Tc? 
(ii) Is superconductivity due to the electron-phonon in- 
teraction or to electron correlation [13? (iii) In the 
case of a phonon mediated pairing which are the rele- 



vant phonon modes [l2t|? (iv) How does the presence of 
electronic donor states (or interlayer states) affect super- 
conductivity |i,[il[ii? 

Two different theoretical explanations has been pro- 
posed for superconductivity in CeCa. In it was noted 
that in most superconducting GICs an interlayer state 
is present at E f and a non-conventional excitonic pair- 
ing mechanismjl^ has been proposed. On the contrary 
Mazin suggested an ordinary electron-phonon pair- 
ing mechanism involving mainly the Ca modes with a 0.4 
isotope exponent for Ca and 0.1 or less for C. However 
this conclusion is not based on calculations of the phonon 
dispersion and of the electron-phonon coupling in CgCa. 
Unfortunately isotope measurements supporting or dis- 
carding these two thesis are not yet available. 

In this work we identify unambiguously the mechanism 
responsible for superconductivity in CeCa. Moreover we 
calculate the phonon dispersion and the electron-phonon 
coupling. We predict the values of the isotope effect ex- 
ponent a for both species. 

We first show that the doping of a graphene layer and 
an electron-phonon mechanism cannot explain the ob- 
served Tc in superconducting GICs. We assume that 
doping acts as a rigid shift of the graphene Fermi level. 
Since the Fermi surface is composed by n electrons, which 
are antisymmetric respect to the graphene layer, the 
out-of-plane phonons do not contribute to the electron- 
phonon coupling A. At weak doping A due to in- plan e 
phonons can be computed using the results of ref. [l5j . 
The band dispersion can be linearized close to the K point 
of the hexagonal structure, and the density of state per 
two-atom graphene unit-cell is iV(0) = /3~^\/87r-\/3-\/A 
with f3 = 14.1 eV and A is the number of electron do- 
nated per unit cell (doping) . Only the E2g modes near F 
and the A'j^ mode near K contribute: 



A = iV(0) 
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where the notation is that of ref. [la| ■ Using this equa- 
tion and typical values of the predicted Tc are 
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FIG. 1: (Color online) (a) Total DOS and DOS projected on selected atomic wavefunctions in CeCa. (b) Band structure of 
CeCa. The size of the dots represents the percentage of Ca component. As a reference the dot at « 0.6 eV at the L-point 
represents 0.95 Ca component, (c) Band structure of *Ca (Red) and C^ (black). The bands have been shifted to compare with 
the CeCa band structure. The directions are given in terms of the rombohedral reciprocal lattice vectors, x is the interception 
between the (rj, —rj, 0) line with the border of the first Brillouin zone {rj = 0.3581). The points L, F, X, T have ri — 0.5. 



order of magnitudes smaller than those observed. As a 
consequence superconductivity in CeCa and in GICs can- 
not be simply interpreted as doping of a graphene layer, 
but it is necessary to consider the GIC's full structure. 

The atomic structure 0| of CaCg involves a stacked 
arrangement of graphene sheets (stacking AAA) with Ca 
atoms occupying interlayer sites above the centers of the 
hexagons (stacking a(3"f). The crystallographic structure 
is R3m [ill where the Ca atoms occupy the la Wyckoff 
position (0,0,0) and the C atoms the 6g positions (x,- 
x,l/2) with x= 1/6. The rombohedral elementary unit 
cell has 7 atoms, lattice parameter 5.17 A and rombo- 
hedral angle 49.55°. The lattice formed by Ca atoms 
in CgCa can be seen as a deformation of that of bulk 
Ca metal. Indeed the fee lattice of the pure Ca can be 
described as a rombohedral lattice with lattice parame- 
ter 3.95 A and angle 60°. Note that the CgCa crystal 
structure is not equivalent to that reported in |ip] which 
has a stacking a/?. In |lOl | the structure determination 
was probably affected by the non-bulk character of the 
samples. 

Density Functional Theory (DFT) calculations are per- 
formed using the PWSCF/espresso code 1 171 and the gen- 
eralized gradient approximation (GGA) [l3| . We use ul- 
trasoft pseudopotentials 19] with valence configurations 
3s^3p^4s^ for Ca and 2s^2p^ for C. The electronic wave- 
functions and the charge density are expanded using 
a 30 and a 300 Ryd cutoff. The dynamical matrices 
and the electron-phonon coupling are calculated using 
Density Functional Perturbation Theory in the linear 
response For the electronic integration in the phonon 
calculation we use a Nk = 6x6x6 uniform k-point 
meshjifiji and and Hermite-Gaussian smearing of 0.1 Ryd. 
For the calculation of the electron-phonon coupling and 
of the electronic density of states (DOS) we use a finer 



Nk = 20 X 20 X 20 mesh. For the A average over the 
phonon momentum q we use a Ng — 4^ q— points mesh. 
The phonon dispersion is obtained by Fourier interpo- 
lation of the dynamical matrices computed on the Nq 
points mesh. 

The DFT band structure is shown in figure ^b). Note 
that the F^X direction and the LF direction are parallel 
and perpendicular to the graphene layers. The K spe- 
cial point of the graphite lattice is refolded at F in this 
structure. For comparison we plot in^c) the band struc- 
ture of CgCa and with Ca atoms removed (Cg*) and the 
structure CeCa with Ce atoms removed (*Ca). The size 
of the red dots in fig. ^b) represents the percentage of 
Ca component in a given band (Lowdin population) . The 
*Ca band has a free electron like dispersion as in fee Ca. 
From the magnitude of the Ca component and from the 
comparison between fig. ^b) and (c) we conclude that 
the CeCa bands can be interpreted as a superposition of 
the *Ca and of the Ce* bands. At the Fermi level, one 
band originates from the free electron like *Ca band and 
disperses in all the directions. The other bands corre- 
spond to the TT bands in Ce* and are weakly dispersive in 
the direction perpendicular to the graphene layers. The 
Ca band has been incorrectly interpreted as an interlayer- 
band not associated to metal orbitals. 

More insight on the electronic states at E/ can be ob- 
tained calculating the electronic DOS. The total DOS, 
fig. n^a), is in agreement with the one of ref. and at 
E/ it is A^(0) = 1.50 states/(eV unit cell). We also report 
in fig. ^a) the atomic-projected density of state using the 
Lowdin populations, pjj{e) = ]^ I]k„ |(</'^|V'kn)P<5(ekn - 

e). In this expression ~ J2riri'^^^^^]v,ri'\4>^') are the 
orthonormalized Lowdin orbitals, \<j>^>) are the atomic 
wavefunctions and S'^_^/ = The Kohn and Sham 

energy bands and wavefunctions are ekn and lipkn)- This 
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FIG. 2; (Color online) (a) and (b) CaCe Phonon dispersion. 
The amount of Ca vibration is indicated by the size of the •, 
of Cz by the size of o, of Cxy by the size of o, of C&xy by the 
size of ▲ and of Caz by the size of T. 



defi-iiition leads to projected DOS which are unambigu- 
ously determined and are independent of the method 
used for the electronic structure calculation. At Ey the 
Ca 4s, Ca 3d, Ca 4p, C 2s, C 2p^ and C 2p^ are 0.124, 
0.368, 0.086, 0.019, 0.003, 0.860 states/(cell eV), respec- 
tively. Most of C DOS at E/ comes from C 2p.n- or- 
bitals. Since the sum of all the projected DOSs is almost 
identical to the total DOS, the electronic states at E/ 
are very well described by a superposition of atomic or- 
bitals. Thus the occurrence of a non-atomic interlayer- 
state, proposed in ref. further excluded. From 

the integral of the projected DOSs we obtain a charge 
transfer of 0.32 electrons (per unit cell) to the Graphite 
layers (A = 0.11). 

The phonon dispersion (wq^) is shown in fig. |21 For a 
given mode v and at a given momentum q, the radii of the 
symbols in fig |5] indicate the square modulus of the dis- 
placement decomposed in Ca and C in-plane (xy, parallel 
to the graphene layer) and out-of-plane (z, perpendicular 
to the graphene layer) contributions. The corresponding 
phonon density of states (PHDOS) are shown in fig. O 
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FIG. 3: (a) Eliashberg function, a F{uj), (continuous line) 
and integrated coupling, X(ijj) (dashed), (b) and (c) PHDOS 
projected on selected vibrations and total PHDOS. 



(b) and (c). The decomposed PHDOS are well separated 
in energy. The graphite modes are weakly dispersing in 
the out-of-plane direction while the Ca modes are three 
dimensional. However the Ca^;;; and the Ca^ vibration 
are well separated contrary to what expected for a per- 
fect fcc-lattice. One CsL^y vibration is an Einstein mode 
being weakly dispersive in all directions. 

The superconducting properties of CgCa can be un- 
derstood calculating the electron-phonon interaction for 
a phonon mode v with momentum q: 



E 

k,n,m 



3kn,k+qm P'^(£kn)^ (Ck+qri 



where the sum is over the Brillouin Zone, 
element is 5k„.k+qm 



(2) 

The matrix 
(knlSV/Suq^lk + qm)/ y^2uj^, 
where Uq^, is the amplitude of the displacement of the 
phonon and V is the Kohn-Sham potential. The electron- 
phonon couphng is A = ^u/Nq = 0.83. We show in 
figEI (a) the Eliashberg function 



1 



2N„^ 



(3) 



and the integral A(w) = 2 J^^duj'a^F{iu')/iu' . Three 
main contributions to A can be identified associated to 
G&xyi Cz and Cxy vibrations. A more precise estimate of 
the different contributions can be obtained noting that 



A = 



[C 



(4) 
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where i, a indexes indicate the displacement in the 
Cartesian direction a of the i^^ atom, [Gqjiaj^ = 
Ek,n,m4s*a5i/3<5(ekn)<5(ek+qm)/[^(0)A^fe], and = 
{lin\SV/Sxqia\li + qm) /^/2- The Cq matrix is the 
Fourier transform of the force constant matrix (the 
derivative of the forces respect to the atomic displace- 
ments) . We decompose A restricting the summation over 
i, a and that over i, /3 on two sets of atoms and Cartesian 
directions. The sets 
resulting A matrix is: 



Ca. 



Cxy Gz G'Axy Cb,z 

/0.12 0.00 0.00 0.00\ 

0.00 0.33 0.04 0.01 (5) 

0.00 0.04 0.27 0.00 

yO.OO 0.01 0.00 0.06/ 

The off-diagonal elements are negligible. The Ca out- 
of-plane and C in-plane contributions are small. For the 
in-plane C displacements, eq. ^ with A = 0.11 gives 
'^c^y.Cxa ~ 0.11. Such a good agreement is probably 
fortuitous given the oversimplified assumptions of the 
model. The main contributions to A come from Ca in- 
plane and C out-of-plane displacements. As we noted 
previously the C out-of-plane vibration do not couple 
with the C tt Fermi surfaces. Thus the coupling to the C 
out-of-plane displacements comes from electrons belong- 
ing to the Ca Fermi surface. Contrary to what expected 
in an fee lattice, the Ca,xy phonon frequencies are smaller 
than the Ga^ ones. This can be explained from the much 
larger A of the Ca in-plane modes. 

The critical superconducting temperature is estimated 
using the McMillan formula 2lj: 



1.2 



exp 



1.04(1 + A) 
A-M*(1-I-0.62A) 



(6) 



where fi* is the screened Coulomb pseudopotential and 
(uj) = 24.7 meV is the phonon frequencies logarithmic 
average. We obtain Tc = IIK, with ^i* = 0.14. We 
calculate the isotope effect by neglecting the dependence 
of ^* on UJ. We calculate the parameter a(X) = 
where X is C or Ca. We get a(Ca) = 0.24 and a(C)''= 
0.26. Our computed a(Ca) is substantially smaller than 
the estimate given in ref. This is due to the fact 

that only « 40% of A comes from the coupling to Ca 
phonon modes and not 85% as stated in ref. |l2| . 

In this work we have shown that superconductivity in 
CeCa is due to an electron-phonon mechanism. The car- 
riers are mostly electrons in the Ca Fermi surface coupled 
with Ca in-plane and C out-of-plane phonons. Coupling 
to both modes is important, as can be easily inferred 
from the calculated isotope exponents a(Ca) = 0.24 and 
a(C) — 0.26. Our results suggest a general mechanism 
for the occurrence of superconductivity in GICs. In or- 
der to stabilize a superconducting state it is necessary to 
have an intercalant Fermi surface since the simple dop- 
ing of the TT bands in graphite does not lead to a sizeable 



electron-phonon coupling. This condition occurs if the 
intercalant band is partially occupied, i. e. when the 
intercalant is not fully ionized. The role played in su- 
perconducting GICs by the intercalant Fermi surface has 
been previously suggested by [i^. More recently a cor- 
relation between the presence of a band, not belonging 
to g raphite, and superconductivity has been observed in 
|l3| . However the attribution of this band to an inter- 
layer state not derived from intercalant atomic orbitals 
is incorrect. 
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